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The number and variability of the neutralizing epitopes targeted by polyclonal 


antibodies in individuals who are SARS-CoV-2 convalescent and vaccinated are key 
determinants of neutralization breadth and the genetic barrier to viral escape’ *. 
Using HIV-1 pseudotypes and plasma selection experiments with vesicular stomatitis 
virus/SARS-CoV-2 chimaeras’, here we show that multiple neutralizing epitopes, 
within and outside the receptor-binding domain, are variably targeted by human 
polyclonal antibodies. Antibody targets coincide with spike sequences that are 
enriched for diversity in natural SARS-CoV-2 populations. By combining plasma- 
selected spike substitutions, we generated synthetic 'polymutant spike protein 
pseudotypes that resisted polyclonal antibody neutralization to a similar degree as 
circulating variants of concern. By aggregating variant of concern-associated and 
antibody-selected spike substitutions into a single polymutant spike protein, we show 
that 20 naturally occurring mutations in the SARS-CoV-2 spike protein are sufficient to 
generate pseudotypes with near-complete resistance to the polyclonal neutralizing 
antibodies generated by individuals who are convalescent or recipients who received 
an mRNA vaccine. However, plasma from individuals who had been infected and 
subsequently received mRNA vaccination neutralized pseudotypes bearing this 
highly resistant SARS-CoV-2 polymutant spike, or diverse sarbecovirus spike proteins. 
Thus, optimally elicited human polyclonal antibodies against SARS-CoV-2 should be 
resilient to substantial future SARS-CoV-2 variation and may confer protection against 
potential future sarbecovirus pandemics. 


Neutralizing antibodies elicited by previous infection or by vaccination 
probably represent a key component of protective immunity against 
SARS-CoV2. Antibodies targeting the receptor-binding domain (RBD) 
ofthe spike protein are thought to dominate the neutralizing activity 
of convalescent or vaccine recipient plasma, and include the most 
potent neutralizing antibodies cloned from individuals who are conva- 
lescent’ ". Nevertheless, additional SARS-CoV-2 neutralizing antibody 
targets include the N-terminal domain (NTD) and the fusion machin- 
ery*®!02, and the full spectrum of epitopes targeted by neutralizing 
antibodies in convalescent or vaccine recipient plasma has not been 
defined. SARS-CoV-2 variants of concern (VOCs) or variants of interest 
(VOIs) encode spike amino acid substitutions" ", some of which confer 
resistance to individual human monoclonal antibodies but have vari- 
able, typically modest, effects on neutralization by polyclonal plasma 
antibodies!9?" ?? Mutated sites in VOCs include those in the RBD, NTD 
and elsewhere, but the numbers and locations of spike substitutions 
required for SARS-CoV-2 to evade the polyclonal neutralizing antibod- 
ies encountered in recipients ofa vaccine or whoare convalescent are 
unknown and are crucial determinants of population immunity. 


Polyclonal neutralizing antibody targets 


Exploiting the fact that SARS-CoV is poorly neutralized by SARS-CoV-2 
convalescent plasma, we compared the sensitivity of HIV-1 pseudotypes 
bearing parental and chimeric spike proteins in which RBD sequences 
were exchanged (SARS-CoV-2(1-RBD) and SARS-CoV(2-RBD); Fig. 1a, 
Extended Data Fig. 1) to neutralization by plasma from 26 individu- 
als from a previously described Rockefeller University COVID-19 con- 
valescent cohort”. The plasma samples were obtained at an average 
of 1.3 months after infection and were selected for high SARS-CoV-2 
neutralization titres (the RU27 plasma panel). Compared with the 
SARS-CoV-2 pseudotype, the SARS-CoV-2(1-RBD) pseudotype was 
less sensitive to neutralization by 21 of the 26 plasmas (median differ- 
ence = 1.8-fold, range 0.5-9.8-fold, P= 0.0005 (Wilcoxon two-tailed 
test); Fig. 1b). Conversely, the SARS-CoV(2-RBD) pseudotype was more 
sensitive than the SARS-CoV pseudotype to all plasmas (median dif- 
ference = 8-fold, range 1.2-75.5-fold, P< 0.0001 (Wilcoxon two-tailed 
test); Fig 1c). Nevertheless, the neutralizing potency of some plasmas 
was hardly affected when the SARS-CoV-2 RBD was replaced by the 
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Fig.1| Neutralizing antibodies in SARS-CoV-2 convalescent plasma targeting 
bothRBD and non-RBD determinants. a, Design of RBD-exchanged chimeric 
spike proteins. b, c, NT;,) for 26 high-titre convalescent plasmas (from the 


SARS-CoV RBD, even though some of those plasmas were minimally 
active against SARS-CoV (for example, RU9, RU1O, RU11 and RUI5; 
Fig. Ib, c). Indeed, plasmas that poorly neutralized SARS-CoV potently 
neutralized chimeric spike pseudotypes with either RBD or the other 
spike domains from SARS-CoV-2 (Fig. 1b, c). For the plasma panel as 
a whole, pseudotype neutralizing potency against SARS-CoV-2 did 
not correlate with potency against SARS-CoV or SARS-CoV-2(1-RBD), 
but did correlate with potency against SARS-CoV(2-RBD) (Extended 
Data Fig. 2a—c). Although altered conformation of the RBD in chi- 
meric spike proteins might affect neutralization, these data indicate 
that while the RBD constitutes a major neutralizing target, substan- 
tial plasma neutralizing activity is also directed against non-RBD 
epitopes. 


Polyclonal antibody-selected variants 


To more precisely map the targets of polyclonal neutralizing antibod- 
ies in individuals who are convalescent, we passaged a recombinant 
vesicular stomatitis virus (rVSV)/SARS-CoV-2 chimeric virus” in the 
presence of each of the RU27 plasmas for up to six passages. rVSV/ 
SARS-CoV-2 mimics the neutralization properties of SARS-CoV-2 (refs. 
+>) but obviates the safety concerns that would accompany such studies 
with authentic SARS-CoV-2. Next-generation sequencing indicated that 
rVSV/SARS-CoV-2 passage in 22 out of the 27 plasmas led to the selection 
of spike mutants (Fig. 2a, Extended Data Fig. 3, Supplementary Table 1). 
For some plasmas, multiple substitutions were selected at distinct but 
proximal sites in viral subpopulations, indicating a dominant neutral- 
izing activity targeting a particular epitope. For other plasmas, substitu- 
tions were enriched in multiple regions of the spike coding sequence, 
suggesting co-dominant neutralizing activities. Six passages of rVSV/ 
SARS-CoV-2 without plasma enriched a small number of substitutions 
that were assumed to represent cell-culture or fitness-enhancing adap- 
tations (for example, T6041), but were distinct from the majority of sub- 
stitutions arising after rVSV/SARS-CoV-2 passage in plasma (Extended 
Data Fig. 3, Supplementary Table 1). Cumulatively, the plasma-selected 
mutations were enriched in specific elements within the NTD, RBD 
and other spike domains (Fig. 2a, Supplementary Table 1). From the 
27 plasma-passaged virus populations, 38 individual mutant viruses 
were isolated by plaque purification; each encoded one, two or three 
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RUI1-27 panel) against pseudotyped HIV-1 virions bearing the indicated spike 
proteins. A median of two independent experiments is plotted. FP, fusion 
peptide; HR, heptad repeat; TM, transmembrane. 


spike substitutions (Fig. 2b) that generally occurred at high frequency 
in the passaged viral populations (Supplementary Table 1). 

We compared the distribution of mutations selected by rVSV/ 
SARS-CoV-2 passage with the RU27 plasma panel with those occurring in 
circulating SARS-CoV-2 populations (Fig. 2a-d).Inboth plasma-selected 
andnaturally occurring sequences, substitutions were enriched in sev- 
eral elements that contributeto the 'supersite' targeted by NTD-binding 
neutralizing antibodies'^P? (Fig. 2a-d). Similar plasma-selected 
and natural sequence variation was also evident in elements 
targeted by class 2 and class 3 RBD-binding neutralizing antibodies”. 
Mutations known to confer resistance to class I RBD antibodies were not 
selected by plasma passage, perhaps reflecting alower than expected 
abundance of class I antibodies in this plasma panel (Fig. 2a-d). Other 
sites, including spike amino acids approximately 680-700 and approxi- 
mately 930, exhibited variation in both plasma-passaged and natural 
variant datasets, but have not yet been demonstrated to be targeted by 
neutralizing antibodies. Nevertheless, the similarity in the distribution 
of natural and plasma-selected sequence variation within the spike 
protein suggests that selection by neutralizing antibodies is a driver 
of divergence in circulating SARS-CoV-2 populations. 

Of the 38 plaque-purified rVSV/SARS-CoV-2 mutants recovered 
following passage in RU27 plasmas, 34 exhibited varying degrees of 
reduced sensitivity to neutralization by the plasma that was used for 
its selection (median 3.1-fold reduced half-maximal neutralizing titre 
(NT.,.), range 0.8-39.3-fold; Extended Data Figs. 4, 5). Nevertheless, for 
37 of the 38 selected rVSV/SARS-CoV-2 mutants, the selecting plasma 
exhibited residual neutralizing activity. We aggregated 13 mutations 
from the plasma-selected viruses based on their effects on plasma 
neutralization sensitivity (Extended Data Figs. 4, 5) and distribution 
throughout the spike protein, generating a single synthetic ‘polymu- 
tant’ spike (PMS) protein sequence, termed PMSI-1 (Extended Data 
Fig. 6a). An rVSV/SARS-CoV-2 derivative encoding these spike muta- 
tions (rVSV/SARS-CoV 2,5.) exhibited resistance to neutralization by 
the RU27 plasma panel that was greater in magnitude and consistency 
than the individual plasma-selected mutants (median 8.0-fold, range 
2.7-52.9-fold; Fig. 2e). Nevertheless, 26 of the 27 RU27 plasmas retained 
residual neutralizing activity against rVSV/SARS-CoV-2pyc,., (Fig. 2e, 
Extended Data Fig. 6b). We conclude that some neutralizing epitopes 
are shared among the convalescent antibodies in high-titre plasmas, 
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Fig.2|Selection of SARS-CoV-2spike mutants by polyclonal antibodies. 

a, Frequencies of amino acid substitutions at each codon of the SARS-CoV-2 
spike protein in two independent rVSV/SARS-CoV-2 populations (1D7 and 2E1), 
determined by Illumina sequencing. Pooled results following selection with 
the RU27 plasma panelare displayed. b, Locations of amino acid substitutions 
in 38 plaque-purified rVSV/SARS-CoV-2 isolates obtained from rVSV/ 
SARS-CoV-2 populations following passage in the RU27 plasmas. c, Frequencies 
of naturally occurring amino acid substitutions (red circles) at each codon of 
the SARS-CoV2 spike protein. The shaded grey bars in a-cindicate shared 


but neutralizing activity against SARS-CoV-2 is clearly polyclonal and 
heterogeneous among individuals with respect to epitope targets. 


Polymutant and variant neutralization 


We generated a panel of HIV-1 pseudotypes bearing the PMS1-1 spike pro- 
tein, asecond PMS protein with a different set of 13 mutations (selected 
based on similar criteria; known as PMSD4), or naturally occurring variants 
or relatives of the SARS-CoV-2 spike protein (Extended Data Fig. 7a, b). The 
panel included several SARS-CoV-2 VOC or VOI spike proteins, and spike 
proteins from sarbecoviruses found in bats (bCoV-RaTG13), pangolins 
(pCoV-GD and pCoV-GX) and previously in humans (SARS-CoV), that 
exhibit varying degrees of sequence divergencefromSARS-CoV-2(refs.”* ?*) 
(see Extended Data Fig. 1 for pseudotype characterization). Todetermine 
sensitivity and resistance to polyclonal SARS-CoV-2 antibodies, we used an 
independentset of 21 randomly selected (Ran1-21) convalescent plasmas, 
and a set of 14 plasmas from recipients who received an mRNA vaccine 
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regions where variation is enriched. d, Comparison of the averaged frequency 
of substitutions observed after passaging rVSV/SARS-CoV-2 with RU27 plasmas 
(centre) and the frequency of sequence changes in natural populations (right), 
projected onto the SARS-CoV-2 spike structure (PDB 6VXX) with positions of 
the RBD and NTD indicated (left). The average frequency of substitutionsina 
15 Aradius is represented using the colour spectrum (scale = 0-20 (centre) 

and 0-9 (right)). e, Neutralization potency of RU27 plasmas against rVSV/ 
SARS-CoV-2 encoding wild type (WT), individually selected mutants or PMSI-1 
spike proteins. A median of two independent determinations is plotted. 


(Vac1-14), inaddition to the RU27 high-titre plasma panel. The PMS spike 
proteins exhibited a degree of neutralization resistance that fell with the 
range of that exhibited by the four SARS-CoV-2 VOCs/VOls and the four 
other sarbecoviruses (Extended Data Figs. 8a—c, 9a-c). Specifically, PMS1-1 
and PMSD4 exhibited neutralization resistance that was greater than 
B.1.1.7 VOC and B.1.526 VOI, similar to P.1 VOC and less than B.1.351.3 VOC. 
PMS1-1 and PMSD4 were more resistant to neutralization than pCoV-GD 
and bCoV-RaTG13, both of which contain a larger number of changes rela- 
tive to SARS-CoV-2 than the PMS spike proteins. Conversely, the pCoV-GX 
and SARS-CoV pseudotypes were more resistant to SARS-CoV-2 conva- 
lescent or vaccine recipient plasma than PMSI-1 and PMSD4 (Extended 
Data Figs. 8a-c, 9a-c). Notably, like PMS1-1 and PMSD4, the B.1.351.3 
VOC that encodes only nine spike mutations relative to SARS-CoV-2 
Wuhan-hu-1 was more resistant to neutralization than sarbecoviruses 
(pCoV-GD and bCoV-RaTG13) that contain a greater number of substitu- 
tions, suggesting that the B.1.351.3 mutations were selected by antibody 
pressure. 
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Fig.3|Neutralization resistance ofpolymutant SARS-CoV-2 spike proteins. 
a, Design of the PMS20 spike protein with 20 antibody-selected and 
VOC-associated mutations. b, Replication of rVSV/SARS-CoV-2 chimaeras 
encoding 2EI (parental) or PMS20 spike proteins in 293T/ACE2cI.22 cells 
infected at a multiplicity of 0.001 and 0.008, respectively. c, Comparative 


Aneutralization-resistant polymutant 


On the basis of the above findings, we attempted to generate a mutant 
SARS-CoV2 spike protein that was minimally divergent from SARS-CoV-2 
Wuhan-hu-1, yet resistant to neutralization by polyclonal convalescent 
and vaccine recipient plasma. Successful derivation of sucha spike pro- 
tein would identify acomplete list of neutralization epitopes recognized 
by polyclonal antibodies. We chose 20 naturally occurring mutations, 
including 8 NTD and 8 RBD changes (Fig. 3a) that either (1) arose in our 
plasma selection experiments (Fig. 2b), (2) occur in VOCs with reduced 
neutralization sensitivity (Extended Data Figs. 7, 8), or (3) arose in our 
previous studies in which human monoclonal antibody resistance was 
selected’”’. Naturally occurring deletion mutations in the NTD (Extended 
Data Fig. 7b), as well as multiple substitutions conferring resistance to 
classes 1, 2 and 3 RBD-binding antibodies'^? were included. An rVSV/ 
SARS-CoV-2 derivative encoding the resulting spike sequence, termed 
PMS20 (Fig. 3a), was replication competent but attenuated compared 
with rVSV/SARS-CoV-2,,,, suggesting that the 20 mutations confer a 
fitness cost (Fig. 3b). Nevertheless, HIV-1 pseudotypes bearing PMS20 
were similarly infectious to those bearing the parental spike protein 
(Extended Data Fig. 1) and were highly resistant to neutralization. Indeed, 
17 of the 21 random convalescent and 8 of the 14 mRNA vaccine recipient 
plasmas produced undetectable neutralization of PMS20 pseudotypes 
(less than 1:50; Fig. 3c). Among the high-titre convalescent RU27 plas- 
mas, 23 of 27 had residual neutralizing activity against PMS20 that was 
reduced by a median of 32-fold compared with the parental pseudo- 


type (range 2.8-114-fold; Extended Data Fig. 9a). We conclude that the 









Polyclonal neutralization breadth 


In contrast to plasmas from individuals who had been infected 
or vaccinated, a panel of plasmas from 14 individuals, termed 
‘infected-then-vaccinated’ (ITV), who had been both infected by 
SARS-CoV-2 and subsequently received mRNA vaccines? retained 
neutralizing activity against HIV-1 pseudotypes bearing the PMS20 
spike (Fig. 4a, b). Indeed, the PMS20 mutations that reduced Ran21 and 
Vacl4 plasma NT;, values by a median of 50-fold (range 5.9-225-fold) 
and 81-fold (range 8.4-229-fold), respectively, caused a median NT; 





Plasma donor 


neutralization potency of randomly selected convalescent (Ran1-21; left) and 
vaccine recipient (Vac1-14; right) plasmas, against Wuhan-hu-1and PMS20 
(b) SARS-CoV-2 HIV-1 pseudotypes. For b andc, the median + range of two 
independent determinations is plotted. 


reduction of only 18.6-fold (range 3.9-100-fold) for the ITV plasma 
panel (Fig. 4a, b). Analysis of chimeric SARS-CoV-2/PMS20 spike pro- 
teins in which the respective RBDs were exchanged (PMS20(2-RBD) 
and SARS-CoV-2(PMS-RBD)) indicated that the relative resistance of 
the PMS20 spike protein to both Ran and ITV plasmas was conferred by 
multiple spike determinants and that the neutralization breadthinthe 
ITV plasmas was due to antibodies directed at both RBD and non-RBD 
determinants (Fig. 4a). In addition to the previously reported potent 
neutralizing activity of ITV plasmas against the B.1.1.7, B.1.525, P.land 
B.1.351.3 VOCs’, the ITV plasmas also potently neutralized B.1.617.2 
(delta), as well as a recently described variant (A.VOI.V2)”* that has 11 
substitutions and 3 deletions in the spike protein, including an exten- 
sively mutated NTD, and is predicted to be resistant to both class 2 and 
class 3 RBD-binding neutralizing antibodies (Extended Data Fig. 10a). 

Plasma from the individuals who are ITV also had substantial neutraliz- 
ing activity against heterologous sarbecovirus HIV-1 pseudotypes, includ- 
ing those that were poorly neutralized by Ran21, Vacl4 and RU27 plasma 
panels and whose RBD and/or NTD sequences are extensively divergent 
from SARS-CoV-2 (Fig. 4c, Extended Data Fig. 10a, b). The median NT, 
values for the ITV plasmas against sarbecovirus pseudotypes were 5,330 
(range 2,369-7,222) for bCoV-RaTG13, 3,617 (range 1,780-6,968) for 
pCoV-GX, 2,605 (range 1,386-3,181) for bCoV-WIV16 and 1,208 (range 
621-2,705) for SARS-CoV (Fig. 4c, Extended Data Fig. 10a). Notably, the 
neutralizing activity of the ITV plasmas against the divergent sarbe- 
coviruses bCoV-WIV16 and SARS-CoV was similar to that found in the 
random convalescent plasmas against SARS-CoV-2 Wuhan-hu-1. Thus, 
the neutralization potency and breadth of polyclonal plasma following 
mRNA vaccination of individuals who had previously been infected with 
SARS-CoV2 seem greater than previously appreciated. 


Discussion 


Theseresults indicatethe presence of abundant neutralizing antibody 
targets on the SARS-CoV-2 spike protein. Our recent analyses further 
suggest that affinity maturation, over months of convalescence, confers 
antibody flexibility and affinity^?? and can impose a requirement for 
multiple viral substitutions for escape from individual neutralizing 
antibodies. Somehuman monoclonal antibodies thus have substantial 
activity against SARS-CoV-2 variants and divergent sarbecoviruses””®. 
Overall, the diversity, maturity and high concentration of neutralizing 
antibodies probably explains why polyclonal plasma from individuals 
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Fig. 4| Neutralization breadth of polyclonal antibodies from ITV 
individuals. a, Comparative neutralization potency (NT;, values) of random 
convalescent (Ran1-2I; left) and ITV (ITV1-14; right) plasmas against HIV-1 
pseudotypes bearing SARS-CoV-2, PMS20 and RBD-exchanged chimeric spike 
proteins. b, Fold differencein NT5,, comparing neutralization of HIV-1 


who have been both infected and subsequently vaccinated could 
effectively neutralize the otherwise highly neutralization-resistant 
PMS20 polymutant, as well as sarbecoviruses that are divergent from 
SARS-CoV-2. While standard mRNA vaccine regimens may be less effec- 
tive than infection at eliciting individual antibody breadth”, it remains 
to be seen whether polyclonal neutralization potency and breadth 
can be achieved using appropriately timed boosting with existing 
SARS-CoV-2 vaccines. If so, existing immunogens may provide robust 
protection against future SARS-CoV-2 variants, and a degree of protec- 
tion against diverse potential future sarbecovirus threats. Conversely, 
PMS proteins encoding numerous neutralization escape mutations 
may represent useful immunogens to broaden the polyclonal antibody 
response elicited by first-generation SARS-CoV-2 vaccines. 
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Methods 


Celllines 

293T (American Type Culture Collection (ATCC) CRL-11268), 293T/ 
ACE2cI.22 and HT1080/ACE2.cl.14 cells? were derived from original stocks 
purchased from the ATCC and were assumed to be authenticated by the 
ATCC. Cells were periodically checked for mycoplasma and retrovirus con- 
tamination by DAPIstaining and reverse transcriptase assays, respectively. 


SARS-CoV-2 and sarbecovirus spike protein pseudotyped 
reporter viruses 

Plasmids pSARS-CoV-2-SA19 and pSARS-CoV-SA19 expressing 
C-terminally truncated SARS-CoV-2 (NC 045512) and SARS-CoV spike 
proteins have been previously described? and were used to construct 
the SARS-CoV-2(1-RBD) and SARS-CoV(2-RBD) expression plasmids 
in which RBD-encoding sequences were reciprocally exchanged. 
A panel of plasmids expressing spike proteins from SARS-CoV-2 VOC 
and VOI was constructed in the context of DSARS-CoV-2-SA19 (R683G)°. 
Substitutions were introduced using synthetic gene fragments (IDT) 
or overlap extension PCR-mediated mutagenesis and Gibson assem- 
bly. All VOCs/VOIs and polymutant spike proteins also included the 
R683G substitution, which disrupts the furin cleavage site and gener- 
ates higher titre virus stocks without significant effects on pseudotyped 
virus neutralization sensitivity (Extended DataFig. 1c, d). The potencies 
with which the plasma neutralized members ofthe mutant pseudotype 
panel were compared with potencies against a ‘wild-type’ SARS-CoV-2 
spike sequence, carrying R683G where appropriate. Plasmids express- 
ing the spike proteins found in the horseshoe bat (Rinolophus affinis) 
coronavirus bCoV-RaTG13 (ref. *’) as well as the pangolin (Manis javan- 
ica) coronaviruses from Guangdong, China (pCoV-GD) and Guanxi, 
China (pCoV-GX)”*” were similarly constructed. Spike sequences 
were codon-modified to maximize homology with the human codon 
usage optimized ofthe pSARS-CoV-2 expressing plasmid VG40589-UT 
(Sinobiological). The19-amino acid truncated coding sequence (CDS) 
of bCoV-RaTG13 (OHR63300), pCoV-GD (CoV EPI ISL 410721) and 
pCoV-GX (CoV EPI ISL 410542) were synthesized by GeneART and 
subcloned into pCR3.1 using Nhel and Xbal and Gibson assembly, and 
referred to as pCR3.1-bCoV-RaTG13-SA19, pCR3.1pCoV-GD-SA19 and 
pCR3.I-pCoV-GX-SAIO9, respectively. Pseudotyped HIV-1 particles were 
generated as previously described". Specifically, virus stocks were 
harvested 48 h after transfection of 293T cells with pHIV-1 GagPol and 
pCCNano/LucGFP (Fig.1) or pNL4-3AEnv-nanoluc (all other figures) 
along with a spike expression plasmid, filtered and stored at -80 °C. 


SARS-CoV-2/sarbecovirus pseudotype neutralization assays 
Plasmas were fivefold serially diluted and then incubated with pseu- 
dotyped HIV-1 reporter virus for 1h at 37 °C. The antibody-pseudo- 
type virus mixture was then added to HT1080/ACE2.cI14 cells. After 
48 h, cells were washed with PBS, lysed with Luciferase Cell Culture 
Lysis reagent (Promega), and Nanoluc luciferase activity in lysates was 
measured using the Nano-Glo Luciferase Assay System (Promega) anda 
Glomax Navigator luminometer (Promega). The relative luminescence 
units were normalized to those derived from cells infected with the 
pseudotyped virus in the absence of plasma. The NT;, was determined 
using four-parameter non-linear regression (least squares regression 
method without weighting) (GraphPad Prism). 


Plasma samples 

Plasma samples were from individuals who were infected with 
SARS-CoV-2 for a mean of 1.3 months before plasma donation’ or from 
individuals who had received MRNA vaccines at various times before 
plasma donation’. A set of 27 plasma samples from individuals infected 
with SARS-CoV-2 with high neutralizing activity who had not been vac- 
cinated’, termed the‘RU27’ plasma panel, were used in VSV-SARS-CoV-2 
selection procedures, while this panel plus a second set of 21 randomly 


selected plasmas (selected at random with blinding to neutralization titre 
or any demographic characteristic) from the same convalescent cohort 
formedthe'Ran2I plasma panel’. A set of 14 plasmas donated by individu- 
als who had received a Pfizer/BioNTech mRNA vaccine formed the ‘Vacl4’ 
plasma panel’. Finally, aset of 15 plasmas from individuals who were con- 
valescent who had received a Pfizer/BioNTech mRNA vaccine between 
6 and 12 months after infection? formed the infected-then-vaccinated 
'ITVI5' plasma panel. The study visits and blood draws were obtained 
with informed consent from all participants under a protocol that was 
reviewed and approved by the Institutional Review Board of the Rock- 
efeller University (IRB no. DRO-1006, ‘Peripheral Blood of Coronavirus 
Survivors to Identify Virus-Neutralizing Antibodies’). 


Selection of antibody-resistant rVSV/SARS-CoV-2 variants 
To select plasma-resistant spike variants, rVSV/SARS-CoV-2/GFP,- 
and rVSV/SARS-CoV-2/GFP,;, were passaged to generate diversity, and 
populations containing 10? plaque-forming units were incubated with 
plasma (diluted 1:50-1:400) for 1h at 37 ?C before inoculation of 2 x 10? 
293T/ACE2cl.22 cells in six-well plates. The following day, the medium 
was replaced with fresh medium containing the same concentrations of 
plasma. Supernatant from the wells containing the highest concentra- 
tions of plasma antibodies that showed evidence of rVSV/SARS-CoV-2/GFP 
replication (large numbers of GFP-positive cells or GFP-positive foci) was 
harvested 24 h later. Thereafter, aliquots (100 ul) of the cleared superna- 
tant from the first passage (p1) were incubated with the same or increased 
concentration of plasma and then used to infect 2 x 10° 293T/ACE2cI.22 
cellsinsix-well plates, as before (p2). In situations where small but expand- 
ing GFP-positive foci were observed, the medium was refreshed at 48 h 
with fresh medium containing no plasma and the virus was harvested 24 
hlater. We repeated this process for up to six passages or until reduced 
neutralization potency for the plasma was obvious, as indicated by visual 
detection of increasing numbers of GFP-positive cells during passage. 
Toisolate individual mutant viruses by limiting dilution, the selected 
rVSV/SARS-CoV-2/GFP,,, andrVSV/SARS-CoV-2/GFP,., populations were 
serially diluted in the absence of plasma and aliquots of each dilution 
added to individual wells of 96-well plates containing 1 x 10^ 293T/ 
ACE2cl.22 cells. Individual viral variants were identified by observing 
single GFP-positive plaques in individual wells at limiting dilutions. 
The plaque-purified viruses were expanded, RNA was extracted and 
spike sequences were determined. 


rVSV/SARS-CoV-2 neutralization assays 

Plasma samples were fivefold serially diluted and then incubated with 
rVSV/SARS-CoV2/GFP,,, and rVSV/SARS-CoV-2/GFP.,, or plaque-purified 
selected variants thereof, for 1h at 37 °C. The antibody-recombinant 
virus mixture was then added to 293T/ACE2.cl22 cells. After 16 h, cells 
were harvested, and infected cells were quantified by flow cytometry. 
The percentage of infected cells was normalized to that derived from cells 
infected with rVSV/SARS-CoV-2 in the absence of plasma. The NT; for 
each plasma was determined using four-parameter non-linear regression 
(leastsquaresregression method without weighting) (GraphPad Prism). 


Sequence analyses 

To identify putative antibody resistance mutations, RNA was isolated 
fromaliquots of supernatant containing selected viral populations or 
individual plaque-purified variants using NucleoSpin 96 Virus Core Kit 
(Macherey-Nagel). The purified RNA was subjected to reverse tran- 
scription using random hexamer primers and the SuperScript VILO 
cDNA Synthesis Kit (Thermo Fisher Scientific). The cDNA was ampli- 
fied using KOD Xtreme Hot Start DNA Polymerase (Millipore Sigma). 
Specifically, a fragment including the coding region of the extracellular 
domain of the spike protein was amplified using primers targeting the 
intergenic region between VSV-M and spike, and the spike intracellular 
domain. The PCR products were purified and sequenced either using 
Sanger sequencing or Illumina sequencing as previously described??. 
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For Illumina sequencing, 1 pl of diluted DNA was used with 0.25 ul Nex- 
tera TDE1 Tagment DNA enzyme (15027865, Illumina) and 1.25 ul TD 
Tagment DNA buffer (15027866, Illumina). Then, the DNA was ligated 
to i5/i7 barcoded primers using the Illumina Nextera XT Index Kit v2 
and KAPA HiFi HotStart ReadyMix (2X; KAPA Biosystems). Next, the 
DNA was purified using AmPure Beads XP (Agencourt), pooled and 
sequenced (paired end) using Illumina MiSeq Nano 300 V2 cycle kits 
(Illumina) at a concentration of 12 pM. 

For analysis of the Illumina sequencing data, adapter sequences 
were removed fromthe raw reads and low-quality reads (Phred quality 
score of less than 20) using BBDuk. Filtered reads were mapped to the 
codon-optimized SARS-CoV-2 S sequence in rVSV/SARS-CoV-2/GFP and 
mutations were annotated using Geneious Prime (version 2020.1.2), 
using a Pvalue cut-off of 10 *. RBD-specific variant frequencies, P values 
and read depth were compiled using Python running pandas (1.0.5), 
numpy (1.18.5) and matplotlib (3.2.2). The parental rVSV/SARS-CoV-2/ 
GFP 2E1and 1D7 sequences each contain two adaptive mutations (F157S 
and R685M for 1D7; D215G and R683G for 2E1), but each was considered 
‘wild type’ for the purposes of the plasma selection experiments and 
were subtracted from the analyses of the sequences. 

Thefrequency of amino acid substitutions during rVSV/SARS-CoV-2 
passage in plasmas was compared with the frequency of global occur- 
rences of changes at each residue on 11 May 2021 (Los Alamos, COVID-19 
Viral Genome Analysis Pipeline; https://cov.lanl.gov/content/index)?.. For 
comparison of SARS-CoV-2 with sarbecoviruses, amino acid sequences 
were aligned with Clustal Omega. Using a python script clone of Simplot 
(https://jonathanrd.com/20-05-02-writing-a-simplot-clone-in-python/), 
the percentidentity relative to SARS-CoV-2 was calculated within a rolling 
window of 100 amino acids, stepping a single residue at a time. 

For three-dimensional sliding window analysis of changes in the spike 
amino acid sequence observed globally and in vitro, the frequency of 
global occurrences of changes at each residue (Los Alamos, COVID-19 
Viral Genome Analysis Pipeline, https://cov.lanl.gov/content/index)*' was 
divided by the average frequency of change at any reside and projected in 
the SARS-CoV-2 spike structure PDB 6VXX” as relative change frequency 
using BioStructMap*’. Alternatively, the averaged frequency of substitu- 
tions observed after passaging rVSV/SARS-CoV-2 with plasma was divided 
by the mean substitution frequency and applied as a three-dimensional 
sliding window over the spike structure. The average frequency of 
substitutions in a 15 A radius is represented using a colour spectrum. 


Reportingsummary 


Further information on research design is available in the Nature 
Research Reporting Summary linked to this paper. 


Data availability 


Data on plasma-selected mutations are provided in Supplementary 
Table 1. Neutralization assay datasets generated during the current 
study are available in the accompanying supplementary source data 
files. Source data are provided with this paper. 
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Extended Data Fig. 1| See next page for caption. 
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Extended Data Fig. 1| Characterization of HIV-1 pseudotypes bearing the 
chimeric, mutant, and variant SARS-CoV-2 and sarbecovirus spike 
proteins. a, Titration of pseudotyped viruses on 293T/ACE2cI.22 cells. 
Chimeric spike pseudotyped viruses inthe upper left panel were built using the 
unaltered SARS-CoV-2A19 and SARS-CoVAI9 spike protein constructs anda 
3-plasmid HIV-1 pseudotyping system (see Methods). The other panels depict 
titration of pseudotypes derived using a furin cleavage site mutant SARS-CoV- 
2A19 spike protein (R683G) and a 2-plasmid HIV-1 pseudotyping system 

(see Methods). b, The same pseudotyped viruses used in (a) were used to infect 


3 different 293T/ACE2 clonal cell lines each expressing a different level of 
ACE2 (MFI= mean fluorescence intensity). c, Titration of pseudotypes bearing 
an unaltered SARS-CoV-2A19 spike protein and a furin cleavage site mutant 
SARS-CoV-2A19 spike protein (R683G) generated using a 2-plasmid HIV-1 
pseudotype system (see Methods). d, Comparative neutralization potency 
(NT; values) of high titer convalescent (RU27) plasmas against HIV-1 
pseudotypes bearing R683G mutant (grey symbols) and unaltered 

(red symbols) SARS-CoV-2A19 spike proteins. For all panels, median + range of 
twoindependent experiments is plotted. 
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Extended DataFig. 2| The RU27 SARS-CoV-2 convalescent plasma panel RU27 plasmas against psudotypes bearing the indicated pairs of spike proteins. 


Simple linear regression was used to calculate R2 and p-values, dashed lines 


contains neutralizing antibodies targeting RBD and non-RBD 
indicate 95% confidence intervals for the regression line. 


determinants. a-c, Correlations of neutralizing potencies of the 
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Extended DataFig.3|Selection pressure on SARS-CoV-2spike exerted by (P2-P6) of two independent rVSV-SARS-CoV-2 populations (1D7 and 2E1), in 
convalescent plasma. Frequencies of amino acid substitutions at each codon each of the RU27 plasmas, determined by NGS sequencing. 
of the SARS-CoV-2 spike protein following the indicated number of passages 
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Extended Data Fig. 4 | Neutralization sensitivity of plasma-selected rVSV/ 
SARS-CoV-2 mutants to RUI-20 plasmas. Infection, relative to 
non-neutralized controls, by plaque purified rVSV/SARS-CoV-2 isolates inthe 
presence of the indicated dilutions of the indicated plasmas from the RU27 


panel. The same plasmas that were used to select the indicated mutants were 
used to determine neutralization potency against the respective plaque 
purified mutants (red) and parental (WT, grey) rVSV/SARS-CoV-2 1D7 or 

2E1 viruses. Median + range of two technical replicates is plotted. 
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Extended Data Fig. 5| Neutralization sensitivity of plasma-selected rVSV/ 
SARS-CoV-2 mutants to RU 21-27 plasmas. Infection, relative to 
non-neutralized controls, by plaque purified rVSV/SARS-CoV-2 isolates in the 
presence of the indicated dilutions of the indicated plasmas from the RU27 


panel. The same plasmas that were used to Select the indicated mutants were 
used to determine neutralization potency against the respective plaque 
purified mutants (red) and parental (WT, grey) rVSV/SARS-CoV-2 1D7 or 

2E1 viruses. Median + range of two technical replicates is plotted. 


K378R 


E484K 
Q493R 
T6041 
F759S 
P792H 


1931M 
V1061M 


o 
N 
© 
= 
Mm 
oO 





RU1 








RU2 









































RU3 RU4 RU5 
b 1 = 1 1 » &— — —1 1 ——$9—9———— 1 —9— —— 
LL. LL. LL LL LL. 
xis 0.1 0.1 = 0.1 = 0.1 a 0.1 a 
Ñ " Ñ Ñ N A 
E E E E E 
5 0.01 5 0.01 5 0.01 5 0.01 5 0.01 
z z z z z 
0.001 0.001 0.001 0.001 0.001 
101 10? 103 104 105 109 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 
Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution 
RU6 RU7 RU8 RU9 RU10 
1 — 1 "—' | 1 1 = a ^ 1 ne —4 
a a a a e a 
LL. LL. LL LL LL 
` 04 fh ` 04 38* 04 ` 04 ` 0.1 pA 
o o o E] "Oo 
o o o o o 
N N N N N 
? ? LE È ? 
S 0.01 5 0.01 5 0.01 5 0.01 5 0.01 
z z z z $ z 
0.001 Tq-9——1—3—1—3 0.001 0.001 T—T—T—T—T—4 0.001 0.001 
101 10? 103 104 105 106 107 101 10? 103 104 105 109 107 101 10? 10? 10^ 105 109 107 101 102 10? 104 105 109 107 101 102 10? 104 105 109 107 
Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution 
RU11 RU12 RU13 RU14 RU15 
1 " 1 1 — 1 v 1 9 
a n. n. . n. a 
LL. LL. LL. LL. LL. 
` 0.1 ` 01 * 0.1 = 0 * 0.1 
ge] "o "o "o "o 
o o o o o 
N N N N N 
E É E E E 
S 0.01 S 0.01 5 0.01 S 0.01 S 0.01 x 
z z z z z ? 
0.001 "T—T—T—T—T—3 0.001 0.001 0.001 0.001 
101 102 103 104 105 106 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 101 102 103 104 105 109 107 
Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution 
RU16 RU17 RU18 RU19 RU20 
1 ————À 1 34 1 NOE ——À 1 oi 1 à 
n. n. n. n. n. ~ 
LL. LL. LL LL LL. 
m 01 m 91 $01 n 94 $001 
N Ñ Ñ N N 
E E Bood * E ? 
5 0.01 5 0.01 s S 0.01 5 0.01 á 5 0.01 
z z z z r z = 
9 
0.001 — DSS S E SN SOE 0.001 0.001 0.001 
10! 10? 10? 104 105 106 107 101 10? 103 104 105 109 107 101 10? 103 10^ 105 109 107 101 102 103 104 105 109 107 10! 10? 10° 10* 10° 10° 107 
Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution 
RU21 RU22 RU23 RU24 RU25 
1 $—— ———À 1 1 —9— — —31 1 «LL 1 
a a a a n 
LL. LL. LL LL LL. 
y f = 04 Ss 04 S 04 a È 04 
ne} o o o o 
o o o o o 
A A A A A e 
Food ! : Fond ; 
E 0.01 E 0.01 E 0.01 5 0.01 5 0.01 
z Z z $ z z 
0.001 0.001 0.001 "T——T—1—1T—T—3 0.001 "T——T——T—1—1— 0.001 
101 10? 10? 104 105 106 107 101 10? 10? 104 105 109 107 101 10? 10? 10^ 105 109 107 101 10? 10? 104 105 109 107 101 10? 10? 104 105 109 107 
Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution Reciprocal plasma dilution 
RU26 RU27 
1 S — —3À 1 ——À 
a a 
LL. LL. 
S $ 
9* 0.1 9* 0.1 
g p rVSV/SARS-CoV-2wr 
T T “© rVSVISARS-COV-2pn5 1-1 
E 0.01 E 0.01 
o o 
z $ z 
0.001 qT—À1—3—3—3À—4 0.001 


1 


0! 10? 10? 104 105 106 107 
Reciprocal plasma dilution 


10! 102 10? 104 10? 109 107 
Reciprocal plasma dilution 


Extended Data Fig. 6 | Neutralization sensitivity of rVSV/SARS-CoV-2 


encoding the PMS1-1 spike. a, Design of the PMS1-1 polymutant spike protein 


with 13 plasma-selected spike mutations aggregated ina single spike. 
b, Infection, relative to non-neutralized controls, by rVSV/SARS-CoV2,ys;, (red) 


plotted. 


and rVSV/SARS-CoV2, (grey) inthe presence onthe indicated dilutions ofthe 
plasmas from the RU27 panel. Median + range of two technical replicates is 
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Extended Data Fig. 7 | Synthetic polymutant and natural variant spike. b, Schematic representation of mutations in naturally occurring 
SARS-CoV-2 spike proteins. a, Design ofthe PMSI-1 and PMSD4 polymutant VOC/VOISARS-CoV-2 spike proteins. 


spike proteins with 13 plasma-selected spike mutations aggregated in each 


Random Convalescent Vaccine recipient Random Convalescent Vaccine recipient 


SARS-CoV-2 SARS-CoV-2 
e PMS1-1 e PMSD4 



















Plasma donor Plasma donor Plasma donor Plasma donor 
b - - 
Random Convalescent Vaccine recipient Random Convalescent Vaccine recipient 
100000 100000 
B.1.1.7 (UK) B.1.526 (NYC) 
10000 10000 
2 os g 
Lr Lr 
Z Z 


pppn 


Ju 






e 
$ 


it 









P.1 (Brazil) 





onor Plasma donor Plasma donor Plasma donor 


RV 
Q 


Plasm 


Random Convalescent Vaccine recipient Random Convalescent Vaccine recipient 
pCOV-GD bCOV-RaTG13 


100000 





SARS-CoV 
10000 
B 
ke 
1000 
$ 
100 
<50 
*—- CO 00 l- SF COO CN SF — 
"v-Cv—-cccc-c 
IINE 
cr x c 
Plasma donor Plasma donor Plasma donor Plasma donor 


Extended Data Fig. 8 | Neutralization potency of random convalescent and against WT (grey symbols) and the indicated SARS-CoV-2 synthetic polymutant 
vaccine recipient plasmas against polymutant, VOC/VOI,andsarbecovirus (a), natural variant (b) or sarbecovirus (c) (red symbol) HIV-1 pseudotypes. 
HIV-1 pseudotypes. a-c, Comparative neutralization potency (NT;, values) of For all panels, median + range of two independent experiments is plotted. 
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Extended Data Fig. 9 | Neutralization potency of high titer convalescent 


HIV-1 pseudotypes. 


a-c, Comparative neutralization potency (NT;, values) of high titer 
convalescent (RU27) plasma against WT (grey symbols) and indicated 


plasma against PMS, VOC/VOI, and sarbecovirus 
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Extended Data Fig. 10 | Neutralization potency of plasma from 
infected-then-vaccinated against VOC/VOI and diverse sarbecovirus HIV-1 
pseudotypes. a, Neutralization potency (NT;, values) of random convalescent 
plasmas (grey symbols) or ITV plasmas (red symbols) against SARS-CoV-2 
prototype or variant or sarbecovirus HIV-1 pseudotypes. Median of two 
independent experiments is plotted. Dashed line indicated median NT,, for 


random convalescent plasmas against Wuhan-Hu-1 SARS-CoV-2. Numbers 
above each scatterplot indicate the median NT;, relative to the medianNT;, 
for Wuhan-Hu-1 SARS-CoV-2. b, Sequence diversity across sarbecovirus spike 
domains; SARS-CoV-2 and the indicated sarbecovirus spike sequences were 
aligned with Clustal and compared using Simplot; the percent identity relative 
to SARS-CoV2 was plotted within a rolling window of 100 amino acids. 
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